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I. 


INTRODUCTION 


The Hne-of-slght component of the solar magnetic field, measured by 
magnetograph In photospherlc lines, has been recorded on an essentially dally 
basis for more than the past quarter century (e.g. data from Mt, Wilson are 
available from 9 August 1959 onwards). Those data have been used as one 
boundary condition for determining potential magnetic field models for the 
region outside the sun (the second needed boundary condition accounts, 
crudely, for the effects of the solar wind by assuming that the magnetic field 
becomes purely radial at a specific spherical surface of radius R in the 
corona, known as the source surface). A comparison of the resulting upward- 
extrapolated magnetic field models with white light photographs of the corona 
has led to the realization that the structure of the corona is, to a large 
extent, controlled by the magnetic field (Newkirk and Altschuler, 1970). 

The primary purpose of the work reported here has been to begin the 
task of determining what the magnetic field patterns are in the region beneath 
the solar photosphere. Part of our Interest in the solar interior has been 
prompted by the recent concentration of attention on "solar seismology" which 
uses the global, free oscillations of the sun to probe Its interior structure, 
Deubner and Gough (1984), For example, that technique now suggests that the 
depth of the solar convection zone is 30% of the solar radius. This is in 
significant contrast both with the range of values of between 15 and 25% of 
the solar radius used in most previous models, as well as the 40% figure * 
favored in some calculations by Gilman (1979). An independent method of 
estimating that depth is clearly of interest, and, at least in the terrestrial 
setting, a magnetic technique Is available (Hide, 1978). It gives a value for 
the radius of the Earth's core-mantle boundary that is In remarkably good 
agreement with the value obtained from terrestrial seismology (Voorhies and 
Benton, 1982). Other topics of interest that can possibly be addressed if 
reliable models of the interior solar magnetic field can be developed include 
elucidating the depth of sunspot magnetic fields and the main electric current 
systems supporting the solar dynamo, plus determining the spin-rate of the 
radiative core, and the strength of the (possibly primordial) magnetic field 
present there. 
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Whereas magnetic probing of the Earth's interior involves the 
relatively straightforward downward extrapolation of a potential magnetic 
field, because the Earth's mantle is believed to be only weakly conducting, 
that physics Is totally Inappropriate for the solar convection zone. Herein 
lies final motivation for this project. The problem of downward extrapolation 
on the sun must rest on a theory of magnetohydrodynamic turbulence, a 
fascinating subject In Its own right, and one In which great theoretical 
strides have been made beginning with the pioneering "mean-field 
electrodynamics" approach of Steenbeck, Radler and Krause (as translated by 
Roberts and Stix, 1971 and developed further by Moffatt, 1978, Parker, 1979, 
Krause and Radler, 1980). Indeed, the following sentence occurs at the end of 
Steenbeck and Krause (1969), c.f. Roberts and Stix, 1971, p. 245: "We would 

like to emphasize once more that, in our opinion, the most Important 
application of our theory lies in the possibility of using it as a probe to 
deduce internal properties of cosmlcal bodies from the magnetic fields 
observed on their surfaces." 

In Section II of this report we sketch our approach to obtaining a 
view of the sub-photospheric magnetic field from magnetic measurements taken 
at the photosphere. The first efforts to implement that theory are then 
described (Section III) and conclusions are summarized in Section IV. 
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II. 


A THEORY FOR THE SUB-PHOTOSPHERIC SOLAR MAGNETIC FIELD 


Suppose that at time t the line-of-slght magnetic field has been 
measured at all points of the solar photosphere, taken to be the sphere 
r 3 Rq and that our goal is to determine the magnetic field (or at least its 
vertical component, for a reason to be supplied shortly) throughout the 
turbulent convection zone. The theory, now developed, is a first step towards 
this difficult goal. 

We first divide the sun into two distinct volumes, the solar 
interior (subscript i) and the solar exterior (subscript e) which are 
separated by the photosphere, r = Rq. It is assumed that electric currents 
flow in r < Rq and r > Rq but that no currents flow in or across the 
photosphere. Then the total magnetic field, $ ( in SI units of tesla, T) at 
the photosphere is obtainable as the gradient of a magnetic scalar potential 
which satisfies Laplace's equation. The solution appropriate to spherical 
coordinates r, 0, <}> ( 0 3 co-latitude, $ 3 east longitude) is a well known 
series of spherical harmonic functions. One obtains the three 
spherical components of t , say 8 r , , and 8^ , by taking the gradient 

of that scalar function. In this way the measured line-of-slght component is 
tantamount to knowing the vertical component of the total magnetic field at 
the photosphere; from here on this last quantity is taken to be the input 
data. 


Next, we decompose the total magnetic field into its parts of 
internal and external origin: 

$ = tL + Z-L , where • (i) 

I 6 

$. = that part of § produced by currents in r < R^ 

§ e = that part of § produced by currents in r > Rq 
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Under the assumption that, classical, low frequency pre-Maxwell electromagnetic 
theory (in which displacement currents are ignored) is appropriate to the sun. 
Ampere's law then implies that 

v x Si a 0 for r > R rt ; v x S a » 0 for r < R* (2) 

Moreover, the absence of monopoles assures us that t is solenoidal and we 
assume that 8^ and & e separately obey the pre-Maxwell equations and appro- 
priate boundary conditions, so that 

V • = 0 , 7 • S 0 a 0 for all r . (3) 

Because is both irrotatlonal and solenoidal in the region on 
and outside the photosphere, its form is the same as that encountered in 
geomagnetism (Chapman and Bartels, 1940) 

°° n R_ n+1 

§i(r,e,<i>,t) = - v(R^ l I (-~) [gJJ(t)cos m4> + hJJ(t)sin m«j>] pj?(e)} (4) 

w n=l m=0 

In this expression, which is valid for r _> Rg , g™(t) , h[^(t) are the time- 
dependent Gauss coefficients, and P^(e) is the Schmidt semi -normalized 
associated Legendre function of order m , degree n . It is next assumed 
that the solar photosphere is sc much closer to strong sub-photospheric 
electric current sources than it is to weak coronal sources, that the Internal 
magnetic field strongly dominates the external magnetic field right at the 
photosphere. Then the measured total vertical magnetic field at r = Rq is 
virtually the vertical component of the internal magnetic field which is 
needed as input to the left hand member of (4). This argument shows that 
standard spherical harmonic analysis can be applied to determine the magnetic 
field of internal origin in the region on and outside the solar surface. 
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Yet, our Interest ts to extrapolate the internal field into the 
solar interior where it need no longer be curl-free. This brings up the 
question of continuity type boundary conditions at the photosphere. 

Imagine the photosphere separates a region of Internal turbulence 
from a region of external qul essence. Then, even though temperature and 
therefore molecular electrical conductivity, a , remain continuous across 
r = Rq , the effective turbulent electrical conductivity changes discon- 
tinuous^ from zero outside the photosphere to some finite, non-zero value 
just inside the convection zone. Across a finite jump in conductivity, there 
can be a jump In vertical electric field associated with a surface electric 
charge, but the tangential components of the electrical field and all 
components of the magnetic field remain continuous. Therefore, with B^ , B Q , 
B^ (actually the components for the Internal magnetic field, but with the 
subscript i usually suppressed from here on) continuous across r = Rq , the 
internal magnetic field can be solenoidal only if 3B r /9r is also continuous 
across r = R because in spherical coordinates equation (3) becomes 


35 r 

ar 


— 8 
r r 


r sin 0 


tw < s e 


sin 0) + 


% 

3<J> J 


(5) 


Therefore, as the photosphere is approached from below, two conditions on the 
internal magnetic field that follow from equation (4) and these two continuity 
conditions are 


at n 

8JR ol . 0, t) = l l (n+1) CgJJ{t) cos m<j> + hJJ(t) sin m<f>3 pJJ(e) (6) 
r 0 n a l m=0 " n n 

3B r (R 0 , 0 , <f>, t)/3r = 

- r' 1 l S (n+1) (n+2) CgJ(t) cos m + hj?{t) sin m*] P>) . (7) 

° n=l m=0 n n n 
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These expressions connect the Internal (vertical) field at the photosphere 
with the Gauss coefficients determined by fitting the line-of-sight measured 
field with spherical harmonics. It should be noted at this point that the 
measurements also determine the tangential components of the internal magnetic 
field as the photosphere Is approached from above, and (by continuity) from 
below; but the vertical derivative fif the tangential components (3B Q /3r and 
3B,/3r) are not continuous across the surface electric charge, so they cannot 

9 

generally be determined In the convection zone without considering that 
distribution Itself. Fortunately, all that Is required for the magnetic 
determination of the depth of the convection zone (as well as some aspects of 
the spin-rate of the radiative core) is the radial component. 

It remains to derive differential equations that determine the 
internal magnetic field In the turbulent solar convection zone. The 
appropriate fundamental equations, written in SI units, with upper case script 
quantities B(T) , E(Vm” 1 ) , JT(Am 2 ) , Q,(Asm 3 ) , t/(ms *) denoting, 
respectively, the (internal) magnetic field, the electric field, electric 
current density, the electric charge density, and fluid velocity relative to a 
sun-centered spherical coordinate system that rotates steadily at the 
Carrington rate, whose period is 27.2753 days, are 


V • £ = 0 , V • I: = £) 

V X s = P 0 ) , V x £ * -3§/3t 


( 8 ) 


5 = o(§ + ft X B) 


Here afSnf 1 ) Is the electrical conductivity, assumed to be a function of 
radius, at most, and e Q » 8.854 x icf i2 As 2 nf 1 T - " 1 , = 4tt • lO -7 T m 1 A 1 

are the dielectric constant and magnetic permeability (taken to be the vacuum 
values) . 
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Next, the total field values are separated Into their large-scale 
mean-field parts (designated by conventional block capitals) and their small- 
scale, turbulent or fluctuating parts (denoted by lower case symbols). Thus, 

'&=$ + £, £»i: + e,5 sl J + l»Q i 3 Q + q»t' a ^ + v, (9) 

where $ a <ti> , £ = <£> , etc. 

and < > denotes the averaging operation. Ideally an "ensemble average", but 
actually a space-time average because of the way In which solar data are 
collected. Quantities linear In the lower case variables have zero mean. 
Introducing this decomposition Into the equations in (8) and then averaging, 
one obtains 


7 4 » & v • £ = cq 1 q 

v x 5 = p 0 3 vx£ = -aS/at do) 

J = a(£ + f 4 + <v x I>>) 

In 0hm { s law, a has been treated as a non-turbulent variable. Subtraction 
of these equations from (0) gives equations for the fluctuating quantltes: 

4* 4 — 1 

7 • b a 0 v • e = Eg q 


7 X 5 = P fl 3 V x e = “3&/8t 


(ID 


j = cr(e + } X 5 + V X I + V X 3 « <{ x S» 


1 1 — 5 



It Is the central task of turbulent dynamo theory to achieve closure 
of these last two sets of equations, for example, by deriving a functional 
relationship between between <v * £> and £ , This is the problem solved so 
persuasively by Steenbeck, Krause and RHdler, We adopt their formalism for 
isotropic mean-field electromagnetism and write, for the 1™ component of 
<\} x (not to be confused with the former use of 1 to designate 
"internal") : 


3B , 

(<V x $>) t = a 1d Bj + U ijk J 


( 12 ) 


where 

Here, 




a 6 


1j 


8 1dk “ 0 e ijk * and 


is a Cartesian coordinate. 


a , 3 are scalar functions of position, determined by the statistics 


of the turbulence. In the more familiar vector notation this gives 


<V X b> « at - sv X '£ , (13) 

Upon substituting (13) into (10) and solving for v x $ (proportional to the 
mean electric current) one obtains 

V X $ = ■ .. Q (£ + ? X 5 + a $) . (14) 

j. + Pq a p 

This turbulent version of Ampere an;! Ohm's law, combined, shows the presence 
of the Important a-effect dynamo term, not present in laminar theory and, more 
Importantly for present purposes, that the effective turbulent diffuslvity is 
not (u 0 a)” 1 as In laminar theory, bqt rather that quantity multiplied by 
1 + lip a 8 , Because £ is the magnetic diffuslvity associated with the 
small scale turbulence, it seems safe to assume that 

« 

0 » (^ a) 1 


(15) 


throughout the convection zone. Consequently, we ignore 1 compared to 
Uq a B in U4) and have, approximately, 

7 X S i 0' 1 (h I X I + a i) . (16) 


This suggests that the appropriate conductivity for the large scale magnetic 
field in the convection zone Is the very small value associated with enhanced 
turbulent magnetic diffuslvity. It would seem that, to a leading approxima- 
tion, large scale electric current does not flow in the convection zone. The 
solar dynamo currents are therefore burled in the small scales. 

This notion can be given mathematical substance by introducing a 
perturbation expansion for the dependence of the mean-field variables on the 
small parameter (1 + a &) -1 that occurs In (14). Let e be a typical 
value of this parameter (or, for a first theory, assume a constant value). 
Then the expansions are 


jt(r, t; a, e) » 5 0 (f, t; a) + e t^(r, t; a) + e 2 $ 2 (r, t; a) + ... 
£(r, t; a, e) = S Q (r, t; a) + e ^{r, t; a) + e Z i^tr, t; a) + ... 


(17) 


Notice that we have not yet made any assumption about the magnitude of a , 
only that the dependence of the mean-magnetic and the electric fields on p 
is analytic so it can be expressed In MacLaurin's series. Upon substitution 
of these expansions into equation (16) and equating the sum of all terms with 
like powers of e to zero, a typical hierarchy of problems arises. The first 
three are 

s°: v • = 0 , v x = 0 (18) 

e*: V • ^ = 0 , V X 8^ = £q + $ x tig + a 8 q (19) 

e 2 : V • § 2 = 0 , V X = ^1 + ^ x $1 + a $1 (20) 
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It Is Interesting that, In the first approximation (for § Q ) we need only 
downward continue the vacuum value of the magnetic field, exactly as If the ^ 

T 

convection zone were an Insulator. The two boundary conditions needed for B^ 

come Immediately from ( 6 ) and (7) with 8 *. replaced by B-. , 

r Up 

All of the numerical work reported In the next section has been 
restricted to this zeroth order approximation. However, on the analytic side 
we have examined the first order correction for finite large-scale current 
flow In the convection zone. For this purpose, It Is convenient to take the 
radial component of the curl of (19) to eliminate in favor of so 

that the first order problem reduces to (for constant 8 but variable a ) 




+ 



1 

r d sin 0 


w < s1 " 9 lr> + 


r z s 1 n z e 3 * 




SB, 


u o a o 


{- 9 ^ + v • [rx (f 


M -FW»oJ 

9 


( 21 ) 


subject to and SB^/ar vanishing at r = R Q . For a specified mean 

motion (e.g., large scale differential rotation) and specified latitudinal 
variation of a effect (e.g., proportional to cos 8 so it changes sign 
across the equator, as required by dynamo theory) and with B Q as determined 
from the zeroth order problem, equation ( 21 ) is a tractable, but complicated 
problem for determining . In future work it is hoped to work out this 
correction. 


He conclude this section by noting that the magnetic method for 
determining the depth of the convection zone is a direct extension of the 
technique introduced by Hide (1978) for the Earth and applied successfully in 
that context by Voorhies and Benton (1982). One calculates the total absolute 
magnetic flux (or pole strength) P crossing the sphere of radius r at 
time t: 
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P(r,t) s /q 71 Jj |B r (r, e, 4>, t)|r 2 sin 0 de d<j> 


( 22 ) 


Note that calculating this integral only requires knowledge of the vertical 
component of the magnetic field. This quantity, proportional to the number of 
field lines crossing the sphere of radius r , has significant variations at 
the solar surface, but, as argued in the proposal leading to this project, the 
very high (molecular) electrical conductivity of the radiative core, together 
with the absence of vertical motion beneath the convection zone, implies that 
P should be independent of time just beneath the bottom of the convection 
zone. In the next section values of P were computed from spherical harmonic 
models of the radial magnetic field (in vacuum approximation) to see if there 
was any tendency for magnetic flux to be frozen at some depth inside the sun. 
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III. 


IMPLEMENTATION OF THE THEORY 


To leading order, the theory described in Section II shows that the 
large scale magnetic field of internal origin remains curl-free within the 
turbulent convection zone. The appropriate technique for downward extrapola- 
tion of the surface observations is then classical spherical harmonic analy- 
sis. That technique has been Implemented in terms of three data sets made 
available to us during this project. Details of those sets are as follows: 

Name of Carrington Time Span of 

Code Data Set Rotation Numbers the Data Set 

A Mt. Wilson (I) 1417-1649 Aug. 1959 - Dec. 1976 

B Mt. Wilson (II) 1620-1717 Sept. 1974 - July 1981 

C Stanford 1642-1742 May 1976 - Nov. 1983 

Data sets A, B, C were supplied, respectively, by Drs. Gordon A. Newkirk, 
Robert Howard, and Todd Hoeksema. Sets A and C had already been fit with 
spherical hamonic expansions truncated at order and degree 9 whereas, set B 
was raw data. 

Initial attention was directed to data set A. The integral in 
equation (22) above was evaluated numerically on a 2° by Z° latitude- 
longitude grid covering the full solar disk using the spherical harmonic 
coefficients of the Mt. Wilson models. The calculations were made for seven 
different spherical surfaces within the convection zone (i"/Rq = 1.00, 0.95, 
0.90, 0.85, 0.80, 0.75, 0.70), and for a random selection of 84 different 
times throughout the 17 year span of the data (but some Carrington rotations 
were selected from each year). The resulting values of P , averaged for all 
of the selected rotations in a given year were then plotted as a function of 
radius ratio, r/RQ . Smooth curves passed through the discrete points of 
such plots would then ideally, cross each other at a distinct radius, thereby 
determining the depth at which absolute magnetic flux becomes invariant with 
time. We would interpret that depth as being close to the bottom of the 
convection zone. 


These first computations produced very few Indications that such a 
depth exists. Instead of two of the curves converging towards each other, as 
depth increases, they mostly diverged, Indicating that when the solar flux 
was, say, relatively small at the surface. It remained small at all depths 
sampled. Actually some of the curves did cross each other, but at angles two 
small to be considered significant. 

In attempting to assess this failure, we came to the belief that 
perhaps the high harmonics In the series were not well determined by the data 
and yet they were dominating the flux integral. Studies of the energy spectra 
for thafse models, as well as the examination of contour plots of the radial 
magnetic field at several different depths confirmed this suspicion. The 
spectra revealed dominance by the eighth or ninth degree terms, much small 
scale structure in the contour maps, and a lot of temporal variability in both 
quantities. In addition, the spurious monopole terms were found to be 
surprisingly large (refer to Figure 1). 

It was, therefore, decided to work with the lower-resolution but 
less noisy data available from the Stanford Solar Observatory, to see whether 
It suffered from the same problems. First, though, the entire span of those 
data (In the form spherical harmonic model coefficients) were passed through a 
low-pass filter (with a cut-off near 2.3 years) designed to remove 
high-frequency noise while retaining the low-frequency solar cycle 
variations. Next, five nearly uniformly spaced Carrington rotations were 
selected (numbers 1648, 1672, 1696, 1716, 1736) and flux computations made, as 
on the Ht. Wilson data. It was again found that there is no depth at which 
absolute magnetic flux appears to be conserved. We believe that the main 
source of difficulty is that noise in the high harmonics is being selectively 
amplified during the downward extrapolation more than that in the low 
harmonics, so that eventually the signal to noise ratio degenerates too much. 
The data should therefore be smoothed .in the spatial domain before downward 
extrapolation. There are sophisticated modern techniques for doing this 
(known as "harmonic splines" and "stochastic inversion" in the geophysical 
literature) but here, we tried only the straightforward technique of 
truncating the spherical harmonic expansions back to lower degress than that 
to which the data were fitted. 
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Figure 1. Normalized energy spectra for 
Carrington rotation 1620. The four 
different traces A, B, C and D refer to 
different depths; 1.0, 0.9, 0.8 and 0.7, 
of a solar radius, respectively. 




The results of this exorcise are encouraging but not definitive, A 
great number of the flux curves now do Indeed cross each other, but the range 
of depths at which this occurs Is substantial. This shows a strong model- 
dependence of the flux conservation depth that will have to be removed before 
we can get results that are really believable. The large amount of numerical 
work on which this rests Is summarized In the three large fold-out matrices. 
Tables III-l, III-2, and II 1-3. The degree of truncation Is N = 5, 7 or 9 . 
Each entry box of each matrix Is at the Intersection of the two Carrington 
rotations being compared (so, close to the main diagonal we are comparing two 
epochs closely separated, while time separation Increases as one moves toward 
the upper right corner). The key on each matrix shows the meaning of some of 
the symbols. The numbers entered Into the boxes show between what radius 
ratios flux was invariant at the two epochs in question. Clearly, as the 
truncation level is lowered from 9 to 7 to 5, there is a substantial growth in 
the number of "successful" events. The region close to the main diagonal 
should be discounted because the pairs of rotations in question are not 
sufficiently distinct. At N = 9, what crossings there are occur at 
unrealistically large depth, but shallower depths, typically down to a few 
tenths of a solar radius, are indicated frequently at N = 7 and N = 5 . 
These results are regarded as encouraging but it would certainly be 
inappropriate to regard them as definitive. There is too much sensitivity to 
the pair of rotations being compared, and probably an overabundance of depths 
too shallow to be believed. 

Our final effort to work with data consisted of a careful fitting of 
spherical harmonics to the second Mt. Wilson data set (Code B). The details 
of this work are included in the Appendix, briefly summarized here. Spherical 
harmonic coefficients, to order and degree 9 (99 internal coefficients exclud- 
ing the monopole) were determined by two different solution techniques (Type I 
and Type II). In solution Type I, all available observations for each 
Carrington rotation were used, with equal weighting, and a single iteration of 
the least squares algorithm was performed. For solution Type II, five itera- 
tions were performed with equal weighting within each iteration, but with each 
data point rejected whose residual of fit at that location was 2.5 times (or 
more) the rms of fit for all data used in the previous iteration. This 
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procedure has produced two different models (with N * 9) for every rotation 
from #1620 to #1717. These models are included In the microfiche packet 
attached to this report. 

To Illustrate the type of behaviour encountered in these data, three 
typical Carrington rotations were selected for detailed study. Rotation #1620 
contains no substantial data gaps, but there are large spatial gradients 
present so the data are noisy (refer to Tables A1-A5 which pertain to solution 
Type I and Tables A16-A20 for solution Type II). Rotation #1625 (see Tables 
A6-A10 for solution Type I and Tables A21-A25 for solution Type II) contained 
data neither excessively noisy nor plagued with data gaps. Rotation #1663 
(Tables A11-A15 for solution Type I and Tables A26-A30 for solution Type II) 
Illustrates a poor distribution of relatively good quality data. The five 
members of these sets of tables reveal, respectively, the distribution of (i) 
the residuals, (1i) their root mean square, (1 1 i ) the number of observations, 
(iv) the llne-of-slght magnetic field calculated from the model, and (v) the 
llne-of-slght magnetic field as observed. 

These calculations will permit culling of the data so as to work 
with a reduced set of the more reliably determined models. We believe that 
such "cleaning procedures" are a necessary step before proceeding further with 
this ’Ine of work. The models for the selected Carrington rotations 1620, 

1625 and 1663 appear in Tables A-31, A-32, and A-33, respectively, for 
solution Type I and in Tables A-34 through A-36 for solution Type II. 

One of the conclusions reached in this study is that, when compared 
with geomagnetic data collected at ground observatories and occasionally by 
satellite, the solar data are far less easy to analyze. A discouraglngly long 
list of potential problems exists and it exacerbates the question of testing 
our theoretical ideas. We do not readily see how to isolate whether the 
theory needs major revision (or refinement) or whether the data are Inadequate 
for the task at hand (or, as is quite possible, both areas require further 
attention). The sources of error which are of concern Include these: 
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1) Instrument error - magnetographs do not provide a perfect 
measure of the 1 Ine-of-sf ght field even within the saturation 
limits, 

2) Instrument saturation - fields In excess of about 200 Gauss fall 
to be delineated. 

3) Absence of strong sunspot field measurements - Insufficient 
light emerges from sunspots to enable spectroscopic 
determination of their strong fields to be recorded, 

4) Absence of limb measurements - line-of-sight fields can weaken 
progressively towards the limb (especially the poles) and fade 
below the limit for detectability. 

5) Data gaps are produced by instrument failure and poor or no 
seeing, 

6) Asynoptic data are gathered but synoptic data are required - 
because reliable data can only be gathered In the region close 
to central meridian, It takes about 27 days to acquire enough 
data for a global model. Time and space aliasing are produced 
by the time evolution during that period. 

7) The tilt of the solar rotation axis by 7* to the ecliptic leaves 
polar caps without data coverage during much of the year, 

8) Differential rotation of the photosphere is not accounted for in 
assigning the location to a given measurement - some locations 
are oversampled while others are not sampled at all, 

9) The source surface , where 6 is assumed to become radial, is 
probably a crude model of solar wind effects because it is taken 
as a sphere of fixed radius. 
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10) Data are, typically, selected off of a grid on which the 
spherical harmonics are not orthogonal , so there Is model 
dependence on truncation level. 

11) There are not trustworthy estimates of modeling and other errors 
available. 



IV, 


CONCLUSIONS 


A concise description of our accomplishments and conclusions are as 

follows: 

A. An approximate method for downward extrapolation of 
line-of-slght magnetic field measurements taken at the solar 
photosphere has been developed. It utilizes the mean-field 
theory of electromagnetism In a form thought to be appropriate 
for the solar convection zone. 

B. A way to test that theory has been proposed. It is based on 
calculating the depth at which the vertical absolute magnetic 
flux remains constant in time and then comparing that value with 
the depth of the solar convection zone as determined by 
helioseismology. Accordingly, existing spherical harmonic 
models of photospheric magnetograph data from Mt. Wilson 
Observatory were used to calculate flux as a function of depth 
for many different solar rotations. 

C. The straightforward application of the lowest-order-theory with 
the complete model fit to these data does not indicate the 
existence of any reasonable depth at which flux conservation is 
achieved. 

D. Spherical harmonic models, fit to data from Stanford Solar 
Observatory at degree N = 9, were first low-pass filtered to 
remove high-frequency noise (ai > 1/2.3 years). Downward 
extrapolation was undertaken on the full models (with N = 9) and 
also on models spatially smoothed by truncating back to N = 7 or 
N = 5. A host of depths at which absolute magnetic flux appears 
to be possibly conserved is then found. However, the range of 
depths is substantial, so there is significant model dependence 
in these results. A definitive frozen-flux depth may well be 
contained in the data, but extracting it in unambiguous fashion 
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will require a considerable amount of Improvement In data 
quality, ?nd perhaps duration of the data span. 

E. Spherical harmonic models, truncated at degree N = 9, have been 
constructed for a span of Mt. Wilson data from Sept. 1974 to 
July 1981 to complement the models previously fit to the Mt. 
Wilson data for the interval Aug. 1959 to Dec. 1976. This gives 
a 22 year continuous record from the same observatory for 
subsequent use. 

F. We speculate that the methods developed and given initial 
testing In this project, may ultimately provide a satisfactory 
magnetic estimate for the depth of the solar convection zone. 
Alternatively, if the value from helioseismology is adopted, 
then the present theory presumably places constraints on the a 
and p turbulence parameters. 
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APPENDIX : SOLAR MAGNETIC FIELD MODELS FROM PHOTOSPHERIC MEASUREMENTS 

The source-free coronal magnetic field may be represented by a 
scalar potential expanded In spherical harmonics 


TT = -vV 


where (see Chapman and Bartels) 
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The expansion may be alternatively written in the form 
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where the coefficients g nm , h nm represent sources internal to the radius a 
and s nm , t represent sources external to a . Then clearly 
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Both Altschuler and Newkirk (1969) and Hoeksema (1984) have utilized scalar 
potential expansions to determine coronal fields from line-of-slght magnetic 
measurements at the photosphere. Because their notations are slightly 
different from each other, we briefly review the mathematical procedure. 


Denote by R^ the radius of the photosphere and R g the radius of 
the source surface In the corona. Then applying the boundary condition that 
the field be radial at R g (V(R g ) = 0) we have from Eq. A.l with a = Rq 
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The scalar potential which describes the coronal field is then 
V = R Q H [(~) n+1 ‘ * n ^°][ A n cos + B n sin m ^ P n (9) 
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In the limit as n gets large, X n approaches zero so that AjJ and bJJ 
approach the internal field coefficients g nm and h nm . The mathematical 
procedure used by Altschuler and Newkirk (1969) and by Newkirk, et al. (1973) 
assumes their limit to be zero. Hoeksema (1984) describes his potential in 
the form 
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so that 

a nm ” ^nm " ^1 - ^nm 
b nm “ ^nm ~ ^1 - x^ B nm 


A. 8 


Hoeksema uses Schmit normalization for his Legendre polynomials and coeffi- 
cients, while Newkirk et al use fully normalized polynomials and coeffi- 
cients. Note that Hoeksema's coefficients are the internal field 
coefficients, while those of Newkirk et al are equal to the internal 
coefficients in the limit as X n approaches zero. 


A- 3 


In thl s appendix we present least squares solar corona field models 
computed using Equation A, 4 and additional photspheric line-of-sight magnetic 
measurements available from the Mt. Wilson magnetograph for Carrington 
rotations 1620 through 1717. These models complement those computed by 
Altschuler, et al (1969) for rotations 1417 through 1649 and by Hoeksema using 
Stanford magneto graph data for rotations 1642 through 1742. The line-of-sight 
measurement at the photosphere Is given by 

B l » B • L 3 cti B r + « 2 B 0 + a 3 


where 


2 2 
= 6 sin 9 cos c - sin 9 + 6 cos 0 

a 2 = 6 sin 0 cos 0 (cos E - 1) - COS 0 


3 6 sin 0 sin e 


A 

and L is the unit line-of-sight vector, e Is the angular Carrington 
longitude increment from the central meridian (CM) to the observation point 
and 
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with a £ the earth-sun distance. Altschuler et al (1969) assume 6 and e 
to be zero. While the measurements are linear functions of the potential 
coefficients, we use an iterative least squares algorithm so that a data 
editing criterion may be applied on sequential iterations. The least squares 
matrix equations 


<Sp n+ i = <A T W A)- 1 A 1 * W Sy n 
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where 


A Is the partial' derivative matrix of the measurements with 
respect to the parameters 

p is the vector of adjusted spherical harmonic potential 
coefficients 

5y Is the vector of residuals, l.e. measured data minus predicted 
value from previous Iteration 

w Is the weight matrix for the measurements 
and the estimate at the (n+l)st iteration Is 


Vl “ p n + %+! 

The weight matrix W Is the Inverse of the noise covariance matrix. We 
assume W to be a diagonal matrix. The matrix (A T W A) 1 is the covariance 
matrix of the estimated coefficients. 

The observation technique for the Mt. Wilson magnetograph and some 
of the limitations of the data are discussed by Altschuler et al (1974). 

Using this data we have computed two different solutions for each Carrington 
rotation from 1620 to 1717; 1) a solution with W set to the unit matrix 
using all data available for that rotation (1080 observations maximum) in a 
single iteration of the least squares equations and 2) a solution computed 
thorugh five Iterations with W set to the unit matrix and with data rejected 
In each iteration step whose residual exceeds 2.5 times the rms sigma of fit 
for the previous iteration (or a minimum of 5 Gauss). 

The solutions for some Carrington rotations suffer from very noisy 
data, while others suffer from poor data distribution. Three representative 
Carrington rotations were chosen; 1620, characterized by noisy data, 1625, 


characterized by good quality data and 1663, characterized by poor data 
distribution, due to gaps In the data. The sigma value and number of data 
points for the three Carrington rotations for solution type I and the fifth 
Iteration of solution type II are given below. 


CARRINGTON 

ROTATION 

SIGMA 

(GAUSS) 

I 

NUMBER OF 
DATA POINTS 

SIGMA 

(GAUSS) 

II 

NUMBER OF 
DATA POINTS 

1620 

3.34 

1056 

0.963 

1026 

1625 

0,817 

1053 

0.799 

1052 

1663 

2.22 

745 

1.192 

722 


For these Carrington rotations spatial distributions averaged on a global 
10° x 10° equiangular grid for the solution residual means, RMS, observed 
values, and computed solution values for the two solution types, together with 
the number of observations within each bin, are presented in Tables A-l 
through A-30. The type II (iterated) solutions clearly remove local regions 
of high Intensity as seen by comparing the rms values in Tables A-2 and A-17. 

The solutions obtained using the first method (which corresponds to 

the procedure used by Altschuler et al (1969)) are presented in Tables A-31 

through A-33. All units are In Gauss. The rms of fit and the number of 

observations are given for each rotation together with the A„ m , B „ coeffl- 

nm nm 

cients and their unit weight standard deviations (from the solution covariance 
matrix) in both Schmidt and fully normalized form. The formal errors for the 
coefficients should be scaled from their unit weight values by the factor 

Mk* rms 

where M is the number of observations and N is the number of solution 
coefficients. The solutions obtained .using the second method which edits data 
so as to exclude anomalous regions of high field strength are presented in 
Tables A-34 through A-36. Solutions (I & II) for all the Carrington rotations 
have been copied onto microfiche and have been included. 




Mean Deviation for solution Type I for Carrington Rotation 1620. Units are 
Gauss with values meaned within a 10° x 10° equiangular grid. 
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RMS Deviation for solution Type I for Carrington Rotation 1620. 
Units are Gauss with values meaned within a 10° x 10° equiangular grid 
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TABLE A-2 



Distribution of observations (multiplied by ten) for Carrington Rotation 1620 
used in solution Type I meaned within a 10° x 10" equiangular grid. 
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TABLE A- 3 


Computed Line-of -Sight Magnetic field measurements from solution 
Type I for Carrington Rotation 1620. Units are Gauss with 
values meaned within a 10° x 10° grid. 
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TABLE A- 4 


Mt. Wilson Line-of-Sight Magnetic field measurements for Carrington 
Rotation 1620 used in solution Type I. Units are Gauss with 
values meaned within a 10° x 10° equiangular grid. 
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TABLE A-5 



Mean Deviation for solution Type I for Carrington Rotation 1625. Units are 
Gauss with values meaned within a 10° x 10° equiangular grid. 
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TABLE A-6 



RMS Deviation for solution Type I for Carrington Rotation 1625. 
Units are Gauss with values meaned within a 10° x 10° equiangular grid. 
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TABLE A-7 


Distribution of observations (multiplied by ten) for Carrington Rotation 1625 
used in solution Type I meaned within a 10" x 10° equiangular grid. 
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TABLE A-8 


Computed Lme-of-Sight Magnetic field measurements from solution 
Type I for Carrington Rotation 1625. Units are Gauss with 
values meaned within a 10° x 10° grid. 
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TABLE A-9 


Mt. Wilson Line-of-Sight Magnetic field measurements for Carrington 
Rotation 1625 used in solution Type I. Units are Gauss with 
values meaned within a 10° x 10° equiangular grid. 
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TABLE A-10 



Mean Deviation for solution Type I for Carrington Rotation 1663. Units are 
Gauss with values meaned within a 10° x 10° equiangular grid. 

0 = 2.22 
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TABLE A-ll 


RMS Deviation for solution Type I for Carrington Rotation 1663. 
Units are Gauss with values meaned within a 10° x 10° equiangular grid 
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TABLE A-12 



Distribution of observations (multiplied by ten) for Carrington Rotation 1663 
used in solution Type I meaned within a 10“ x 10° equiangular grid. 
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TABLE A-13 



Computed Line-of-Sight Magnetic field measurements from solution 
Type I for Carrington Rotation 1663. Units are Gauss with 
values meaned within a 10° x 10° grid. 
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TABLE A-14 



Mt. Wilson Line-of-Sight Magnetic field measurements for Carrington 
Rotation 1663 used in solution Type I. Units are Gauss with 
values iteaned within a 10° x 10° equiangular grid. 
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TABLE A-15 



Mean Deviation for the fifth iteration of solution Type II for 
Carrington Rotation 1620. Units are Gauss with values meaned 
within a 10° x 10° equiangular grid. 
a = .963 
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TABLE A-16 


RMS Deviation for the fifth iteration of solution Type II for 
Carrington Rotation 1620. Units are Gauss with values meaned 
within a 10° x 10° equiangular grid. 
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TABLE A-17 


Distribution of observations (multiplied by ten) for Carrington 
Rotation 1620 used in the fifth iteration of solution Type II 
meaned within a 10° x 10° equiangular grid. 
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TABLE A-18 



Computed Line-of-Sight Magnetic field measurements from the fifth 
iteration of solution Type II for Carrington Rotation 1620. Units 
are Gauss with values meaned within a 10° x 10° grid. 
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TABLE A-19 



Mt. Wilson Line-of-Sight Magnetic field measurements for Carrington 
Rotation 1620 used in the fifth iteration of solution Type II. Units 
are Gauss with values meaned within a 10° x 10° equiangular grid. 
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TABLE A-20 



Mean Deviation for the fifth iteration of solution Type II for 
Carrington Rotation 1625. Units are Gauss with values meaned 
within a 10° x 10° equiangular grid, 
a = .799 
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TABLE A-21 


RMS Deviation for the fifth iteration of solution Type II for 
Carrington Rotation 1625. Units are Gauss with values meaned 
within a 10“ x 10° equiangular grid. 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

Qy 


OO^OChHlO^OOOfOOQOOOIOOvOOOOOOOOOQNOOtVChtfOO 

fH rH NNN K) 


o^iriin^'J , K)OK)NiHOQNC 4 <MC'C , K 7 f^fH^fHiriiHM'a-CJro»OfMWK)ON»HCvi<r 
f>% fH HH 


QL 0 ' 3 * 0 >mUl(MOLns 0 K)in^ 0 fHS 3 -fHKl'CrC\N‘ 0 ^’f 0 ^)C\^ 0 a^‘a'sr{MMvJ , «?C 0 M ?0 
VO fH 


Or^CSlK)OK)COS 3 , 0 »nmO<SlOOU 1 r^CJiHf>tf>^MCONFH«IOfHNin'a'K 7 ^'ONrH 
lf> rH rH fH 


>r^*a - covTfor^r>iOK>rN.*»TsrfoinsDrHirjMsTfOio , a'<Njror^^r>ofo^uih-^>*HK>f^vi) 


QK)rHO^O' 3 ’lOfO'Oin*J , NM^OOCONOlOlOfVlOfO' 00 >NLnO'NTf>.CO*a'f^<rvr^lO 
K)H H H 


apjin«jiOi/^iovrMh.iHr^iHoofHarO'Oio^ro>K»fs.^'TNvoofs.rH\o>o^h-oK> 

CM iH H »H CM rH 


OC 0 'CrOC>fHS 3 ’ir|K)*JT^ k 0 C 0 >TK>irur|CC f, '<Hf^s 3 'M?'T' 0 VrsTf 0 l 0 fHiHC>CQ'^iHCM ‘0 
rH fH rH HrtNNHHNN rH 


jvj*r^<M»Ain»o>oin^O'oooLncMco(NJ^<McCM^Dfo^fs.minco^oeof>.ovcooMjineoift 


Of0OCMt0f^C7'MrOa'0tn'0^CVJCr»\0<ri0r0f0CMCM'J , r0C0CMsTC0Oc0fHr^f0CM>0^3 
rH rH rH rH CM fO rH rH CM CM K) K) 


ovT^srNsT'ofOhO^ocnrN.NvovTfftvT*'>inv 0 »Hr^s 3 'N , >vr’a , r , w©c | O(M^or , '-'< 3 'o> l co'o 

CM 


OSTK>K)rH' 0 inCVJrN<Min^^h.C 0 'a'' 3 , O<Min 00 rHCMKlK>»HtninC<Mf 0 omCMCMl 0 «M 
K> rH rH 

I 

OOMK )'0 0 'NH' 3 , OmHMKJSHOf 0 lfl<f « 0 r 0 f 0 lf|O' 0 f 0 OHHHartOOHC' 

rH rH rH rH fH rH 


omoofOvooinoo'a'rHin\or»Jfors.io\orHK)toN , >cMroro^*oincoK>rH<MNominsrKiw 

lo rH rH 


©COO*4*STOOrHfOOrHCMOOrHO*a'CMO'iCMOOOWlvr*OCC(MiHr'-rHlOO>COrHir|OCMn 
vO H rH rH iH 


OHinHOsM , 'OOv^O(\jHHOHHOWK)STKlN«OinO'OHO«(SlONH)H(VJwJO' 
f*. rH CM rH rH <M rH rH (M rH rH rH 


Doaooooooooooooooaoooooooooooooooaaoo 

CO 


rH(M> 0 »crm^ 0 f^C 00 N©rH<M> 0 'J*mM>r *00 0 '©rHCMr 0 vTin^ 0 r^ 00 CNOrHtMh 0 'Tin ^0 

rHHrHrHrHrH»HrHrHrHCMCMNCMCMCMCMCMCMNfOK7fOK)K)K)K1 


A-28 


TABLE A- 22 


Distribution of observations (multiplied by ten) for Carrington 
Rotation 1625 used in the fifth iteration of solution Type II 
meaned within a 10° x 10° equiangular grid. 
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TABLE A-23 



Computed Line-of-Sight Magnetic field measurements from the fifth 
iteration of solution Type II for Carrington Rotation 1625. Units 
are Gauss with values meaned within a 10° x 10° grid. 
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Mt. Wilson Line-of-Signt Magnetic field measurements for Carrington 
Rotation 1625 used in the fifth iteration of solution Type II. Units 
are Gauss with values meaned within a 10° x 10° equiangular grid. 
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Mean Deviation for the fifth iteration of solution Type II for 
Carrington Rotation 1663. Units are Gauss with values meaned 
within a 10° x 10° equiangular grid, 
a = 1.192 
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RMS Deviation for the fifth iteration of solution Type II for 
Carrington Rotation 1663. Units are Gauss with values meaned 
within a 10° x 10° equiangular grid. 
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Distribution of observations (multiplied by ten) for Carrington 
Rotation 1663 used in the fifth iteration of solution Type II 
meaned within a 10° x 10° equiangular grid. 
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Computed Line-of-Sight Magnetic field measurements from the fifth 
iteration of solution Type II for Carrington Rotation 1663. Units 
are Gauss with values meaned within a 10° x 10" grid. 
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